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ABSTRACT 

We present 4 model series of the CODEX dynamical opacity-sampling models of Mira variables 
with solar abundances, designed to have parameters similar to o Cet, R Leo and R Cas. We 
demonstrate that the CODEX models provide a clear physical basis for the molecular shell 
scenario used to explain interferometric observations of Mira variables. We show that these 
models generally provide a good match to photometry and interferometry at wavelengths 
between the near-infrared and the radio, and make the model outputs publicly available. These 
model also demonstrate that, in order to match visible and infrared observations, the Fe-poor 
silicate grains that form within 3 continuum radii must have small grain radii and therefore 
can not drive the winds from 0-rich Mira variables. 

Key words: stars: AGB and post-AGB, stars: atmospheres, stars: variables: other, stars: mass 
loss 
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1 INTRODUCTION 

] Asymptotic Giant Branch (AGB) stars represent the final fusion- 
: powered stage in the evolution of solar-type stars, and the engine 
by which the vast majority of the material in o ur Galaxy i s recy- 
cled from stars back to the interstellar medium (Gail 2003). Mira 
variables represent the final stage in AGB evolution before they be- 
come dust-enshrouded and difficult to observe. They are so bright 
in the infrared that they can be used as competitive extragalactic 
distance indicators and pro bes of star formation history (Rejkuba 
I2OO4I : iMenzies et al.ll2010)) . They are also unique amongst stellar 
classes in their opportunity for detailed observations: light curves 
that differ in shape and amplitude in different bandpasses, pho- 
tospheres that show different sizes and structure as a function of 
wavelength, velocity-resolved motions and complex spectra. 

The observational literature on Mira variables is very ex- 
tensive, as partly detailed below in Section |4] In order to make 
sense of these observations, comprehensive models are required 
that link physical parameters to pulsation, and pulsation to ob- 
served p roperties and mass- l oss. The previous g eneration of mod- 
els (e.g. iHofner et al.lll998l : iHofmann et al.|[T998h suffered from 
grey or mean-opacity like approximations in their radiative trans- 
fer codes, so were not ideally suite to interpreting many observed 
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properties, for example visible brightness or high-res olution spec- 
trosco py. Th e next generation of models (Upsalla lHofner et al.l 
l200l cODEx llreland et"al]|2008h are better suited for modern ob- 
servational comparisons, but extensive grids have not yet been pro- 
duced, both because of solvable but difficult computational issues 
and the lack of a clear calibration for model parameters. 

Here we present four physical model series for M-type Mira 
variables, as a first step in tuning and testing the CODEX models 
to derive physical parameters of Mira variables from observations, 
and to gain physical insight into the dominant physical processes 
in Mira-variable atmospheres. Observational predictions including 
predictions for infrared interferometry of the model series are made 
available online so that new observations can easily be compared to 
these models. 

One of these model series has parameters based on o Cet, one 
based on R Leo and two are based on R Cas with different as- 
sumptions. For all models, only the pulsation period is guaranteed 
to match observations, and in this paper we aim to examine the 
other model outputs in order to more closely target model series to 
real stars in future papers. The model construction was described 
in Paper I (Ireland et al. 2008): they begin with input parameters 
of mass, luminosity and coinposition, with three other free param- 
eters being microturbulent velocity, mixing length (Om) and turbu- 
lent viscosity (a^). Pulsation is self-excited (i.e. it occurs sponta- 
neously), and the temperature of all layers is calculated by solving 
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the conservation of energy equation via a detailed opacity-sampling 
metliod. 

Details of the parameter choices for the model series are given 
in Section[2l as well as basic comparison of the model light curves 
to observations. Model predictions for spectra, in particular the ef- 
fect of extension on spectra, are tested in Section[3] and the models 
are compared to observations of o Cet including infrared and radio 
interferometry in Section |4] In Section |6] we compare the model 
structures to previously published ad-hoc molecular shell papers, 
and in Section|7]we discuss the mass loss rates of the models and 
the driving mechanisms. In Section|5]we discuss the effect of input 
parameters on the models, and the possibility for better calibrat- 
ing the input parameters so that, e.g. the mass of individual Miras 
could be inferred from models. Finally, in Section [8] we conclude 
and discuss plans for future work. 



2 MODEL PARAMETERS AND DESCRIPTION 

A detailed description of the model construction was given in 
Paper I. Briefly, the models consist of self-excited grey models 
that determine the atmospheric pressure stratification and lumi- 
nosity. The temperature profile is then re-iterated using an opac- 
ity sampling code with 4300 wavelength points, assuming radia- 
tive and local thermody namic equilibrium. Dust formation follows 
llreland & Schola yOOg), except that we drop the Rayleigh approx- 
imation, instead replacing it by a smooth fit to the Mie approxima- 
tion of spherical grains, weighting the Rayleigh scattering by: 



(JM{a,\) 



on{a,\){l + A.^{^fr 



(1) 



where a is the grain radius, A the wavelength of radiation, 
a the scattering cross-section and or the scattering cross-section 
in the Rayleigh approximation. This cross-section om is the total 
cross-section weighted by (1-^), where ^l is the impact parame- 
ter. This weighting ensures that the radiative acceleration on dust is 
correct in the presence of forward-scattering. 

The choice of free parameters was only briefly discussed in 
Paper I in the context of the o 5 4 model series, based on the param- 
eters of the prototype Mira variable o Cet. The input parameters 
for all 4 model series presented here are given in Table [T] and the 
reasons for their choice give in sub-sections below. 

The behavior over lO** days of each nonlinear pulsation model 
series is shown in Figures 1-4. As each model series runs for many 
cycles, we chose only a few typical cycles for detailed examination. 
In each of these cycles, ~I0 representative models were extracted 
and their velocity and pressure structures used as input to the model 
atmosphere code which, after temperature iteration, give spectra 
and centre-to-limb variation (CLV) for the models. The cycles dur- 
ing which models were extracted are shown as shaded regions in 
Figures \-4 and the actual models extracted are shown as circles in 
the top panels of these figures. The instantaneous physical param- 
eters and shock front locations for the chosen models are given in 
Tables [3}{8] (available in their entirety in the on-line version of the 
journal). We note that the model radius is defined in these tables 
as the radius where the Rosseland optical depth is unity, and the 
effective temperature is defined by this radius. 

For all 4 m odel series, we assume solar element abundances 
from lOrevesse e t al. ( 1996). This is near the mean abundance ob- 
served for stars in the s olar vicinity with ages of 3-6x10^ years 
JEdvardsson et alj|l993h . Our red giant model with mass 1.35 M© 
has an age of ~3 x l(f years and the red giant model with mass 
I.I Mq has an age of ~6xl0^ years according the isochrones 



Table 1. Parameters of 4 model series. The mass M, luminosity L. metal- 
licity Z. mixing-length parameter cim and turbulent viscosity parameter 
Op are input pai'ameters, and the parent-stai' radius Rp and period P are 
derived parameters. 



Name M 



Rv 



Puno 



o54 


1.1 


5400 


0.02 


3.5 


0.25 


216 


330 


261 


R52 


1.1 


5200 


0.02 


3.5 


0.25 


209 


307 


243 


C50 


1.35 


5050 


0.02 


2.0 


0.24 


291 


427 


408 


C81 


1.35 


8160 


0.02 


3.5 


0.32 


278 


430 


374 



M=l.l M 



L=5400Lg Z = 0.02 P = 330 days 




2000 



4000 6000 

time(days) 



8000 



Figure 1. The luminosity (top panel), effective temperature (middle panel), 
and the radii of a representative selection of mass zones (bottom panel) plot- 
ted against time for model o 5 4 . The red line in the bottom panel shows the 
position of the point where the grey approximation optical depth Tg = ^. 
Tcff in the middle panel is here defined as the temperature where Tg = 2/3, 
which is close to the effective temperature oc (L/i?,^)^/* of the non-grey 
atmospheric stratification. The shaded regions show the time intervals in 
which models were selected for detailed atmospheric model computation. 
The selected models are circled in the top panel. The mass, luminosity, 
metallicity and period of the nonlinear pulsation model are shown at the 
top of the plot. 



of lOirardi et alJ ( l2OO0 h. However, as shown by Isdvardsson et al.l 
( 119931) , there is a scatter in [Fe/H] from about -0.5 to -1-0.3 for stars 
of this age, so that the Mira stars we are aiming to model could 
have somewhat different abundances to those we have adopted. 

The adopted value for the mixing length in units of pressure- 
scale height for the o54, R52 and C81 series, am. = 3.5, is unusu- 
ally high compared with models of more compact stars. However, 
this is not unreasonable, as detailed hydrodynamic calculations of 
stellar convection often su ggest values for am in th e range 1-4 (see 
the parameter summary in lMeakin & Ametll2007h . Decreasing the 
mixing length in models makes heat transport more difficult in the 
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Figure 2. The luminosity and mass-zone positions in the R52 series, show- 
ing the locations of the phases chosen for detailed model computation. 
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Figure 3. The luminosity and mass-zone positions in the C 8 1 series, show- 
ing the locations of the phases chosen for detailed model computation. 



outer convective layers, causing the model star to expand. There- 
fore mixing length has been used by us as a way to produce the 
correct period, given a luminosity L. However, this procedure is 
always ambiguous for a given field Mira, because the distance and 
hence M is always uncertain by ~10% or more. 

The value of q^ can be adjusted to give the correct pulsation 
amplitude. Alternatively, within limits, M can adjusted to give the 
correct pulsation amplitude. Thus far in our models series, we have 
assumed values of M based on e.g. typical masses as a function of 
period and then used a^ to tune the pulsation amplitude. By com- 
paring these model series with observations, we aim to develop a 
preferred values for Qm and a^, or at least a preferred prescrip- 
tion for choosing q,„ and q„ as a function of the physical input 
parameters M, L (a proxy for evolutionary state along the AGB) 
andZ. 

Although our models have been generally based on the stars 
o Get, R Leo and R Cas, there are a range of parameters that are 
consistent with these Miras, and a detailed comparison with ob- 
servations will inevitably reveal where the differences lie. In the 
following section, we aim to discuss the chosen parameters for 
each model series, compare predicted light curves to observed light 
curves and draw preliminary conclusions as to whether changes to 
physical input parameters could improve model fits. We choose the 
V, J and K band for comparison to observations: V band because 
of the wealth of observational data, and J and K because the model 
predictions are most reliable in these band-passes. For this com- 
parison, we add 0.5 magnitudes to the predicted V-band fluxes, to 
account for non-LTE effects as computed in the models at selected 
phases in Paper I, where the correction needed varies between 0.3 
and 1.0 magnitudes. 



Table 2. Parameters of the extended o 5 4 cycle, including the position of 
the shock fronts. 



Model 


Phase 


L 


R 


Tcft 


Si 


S2 


S3 






(i©) 


(Rp) 


(K) 


(Rp) 


(Rp) 


(Rp) 


260820 


-0.20 


4050 


0.90 


3299 


2.82 


0.91 




260960 


-0.10 


4240 


0.85 


3434 


2.85 


0.87 




261140 


0.00 


8013 


0.98 


3761 


2.86 


0.99 




261320 


0.10 


8102 


1.20 


3398 


2.87 


1.28 




261460 


0.19 


7420 


1.34 


3154 


2.83 


1.50 




261620 


0.31 


5830 


1.41 


2898 


2.74 


1.70 




261740 


0.40 


4462 


1.42 


2697 


2.64 


1.82 




261860 


0.50 


3640 


1.34 


2640 


2.47 


1.91 




261940 


0.60 


2440 


1.40 


2333 


2.31 


1.93 




262160 


0.70 


3450 


0.97 


3055 


2.07 


1.89 


0.97 


262360 


0.80 


4110 


0.90 


3325 


=> 


1.82 


0.90 


262600 


0.90 


4355 


0.85 


3462 




1.74 


0.86 


263160 


1.00 


8428 


0.99 


3786 




1.61 


1.01 


263740 


1.10 


8420 


1.20 


3439 




=> 


1.37 



2.1 o 5 4 5400 Lq series 

The 5400 Lq model series for o Get (P = 332 days) has model 
parameters chosen to match the luminosity obtained from the J 
and K photometry of Iwhitelock et alj J200P ) and the lKnapp et al.l 
(2003) revised HIPPARCOS parallax. The mass of 1.1 M© was 
adopted since an analysis of the G alactic scale height of Mira vari- 
ables bv ljura & KleinmanrJ ( 119920 suggest a mass of '--^l.l M© for 
Mira variables with periods from 300-400 days, while a study of 
the population of Miras by Wyatt and Gahn (1983) gives a progen- 
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Table 4. Parameters of the 4-cycle continuous phase coverage o54 cycle, including the position of the shock fronts. [To be available on-line only] 



Model 


Phase 


L 


R 


Toff 


SI 


S2 


S3 


S4 


S5 S6 






(Lq) 


(Rp) 


(K) 


(Rp) 


{Rp) 


(Rp) 


(Rp) 


(Rp) (Rp) 


285180 


-0.20 


3979 


0.92 


3263 


3.54 


0.92 








285320 


-0.11 


4107 


0.86 


3397 


4.06 


0.87 








285500 


0.00 


7288 


0.93 


3769 


4.21 


0.94 








285700 


0.10 


8236 


1.17 


3468 


4.36 


1.22 








285860 


0.20 


7623 


1.33 


3189 


4.52 


1.47 








285980 


0.31 


6270 


1.40 


2952 




1.70 








286060 


0.41 


4513 


1.43 


2692 




1.87 








286100 


0.49 


3845 


1.40 


2614 




1.95 








286140 


0.59 


2548 


1.53 


2258 




2.02 








286320 


0.70 


3204 


1.00 


2959 




2.00 


1.00 






286520 


0.80 


4092 


0.91 


3291 




1.93 


0.92 






286700 


0.90 


4218 


0.86 


3415 




1.83 


0.87 






287000 


1.00 


7188 


0.93 


3755 




1.64 


0.94 






287560 


1.10 


8431 


1.17 


3487 




1.38 


1.25 






287740 


1.20 


7734 


1.32 


3204 




=► 


1.62 






287820 


1.30 


6304 


1.40 


2963 






1.95 






287880 


1.40 


4658 


1.43 


2720 






2.18 






287940 


1.51 


3820 


1.38 


2629 






2.38 






287980 


1.61 


2661 


1.47 


2326 






2.52 


1.13 




288140 


1.70 


3342 


1.00 


2984 






2.61 


1.00 




288320 


1.80 


4034 


0.91 


3279 






2.67 


0.92 




288460 


1.90 


4159 


0.86 


3407 






2.71 


0.87 




288620 


2.00 


7340 


0.93 


3770 






2.71 


1.37 




288820 


2.10 


8170 


1.17 


3461 






2.71 


1.22 




289020 


2.20 


7495 


1.32 


3181 






2.69 


1.46 




289240 


2.30 


6174 


1.40 


2948 






2.66 


1.66 




289440 


2.40 


4655 


1.42 


2730 






2.63 


1.80 




289620 


2.49 


3918 


1.36 


2671 






2.46 


1.91 




289740 


2.59 


2525 


1.42 


2336 






2.26 


1.97 




289920 


2.70 


3064 


0.99 


2932 






=► 


1.99 


0.99 


290120 


2.80 


4073 


0.91 


3290 








2.01 


0.92 


290360 


2.90 


4301 


0.86 


3434 








2.08 


0.87 


290740 


3.00 


7795 


0.95 


3796 








1.97 


0.96 


291500 


3.10 


8412 


1.17 


3474 








2.06 


1.29 


291740 


3.19 


7759 


1.32 


3210 








=^ 


1.65 


291800 


3.31 


6097 


1.40 


2935 










2.19 


291820 


3.41 


4451 


1.43 


2688 










2.55 


291840 


3.55 


3438 


1.40 


2544 










3.07 


291860 


3.61 


2813 


1.35 


2469 










3.24 1.13 



itor mass estimate of 1.18 M© for o Cat. We give this model series 
the designation o54. 

Figure [5] shows the light curves derived from 4 cycles of the 
o54 model series between times 8000 and 10* days in Figure[T] 
To fit the K-band maximum flux, o Cet needs to be at a distance 
of 107 pc . This corresponds to a parallax of 9.3 mas, consistent 
within 2(7 of the late st HIPPARCOS value of 10.91±1.22mas from 
Ivan LeeuwenI d2007r) . 

The models fit the general light-curve shape and amplitude, 
but they are ~0. 1 magnitudes too blue in J-K colour, too blue in 
V-K by ~1 magnitudes (hence too bright in V near maximum by 
~1 magnitudes), and these 4 cycles do not reproduce the cycle-to- 
cycle scatter in J and K magnitudes near minimum. 

We also computed models for more compact atmosphere cy- 
cles (between times 1000 and 2000 days in Figure 1) and more ex- 
tended atmosphere cycles (between times 3500 and 4300 days), but 
J, H and V fluxes were not noticeably different. The atmospheric 
structures for the various models are shown in Figure [T] The con- 
sistently blue J-K and V-K colours suggest that Teff is too high 
in the model (see Section 5 for more discussion). 

Although the basic model properties (L, R, Tcs) and the near- 



continuum fluxes are similar from cycle to cycle, high-layer obser- 
vational features differ depending on the upper atmosphere struc- 
ture. The reason for this can be seen in Tables [2] through [4] where 
upper layer shock front positions are not repeatable from cycle to 
cycle. For example, in Table|2] at phase -0.2, the upper strong shock 
from a previous cycle is at 2.82 Rp. Exactly 1 cycle later, at phase 
0.8, the upper shock has just merged with a lower shock at 1.82 Rp. 

2.2 R5 2 5200 Lq series 

The period of R Leo is slightly shorter than that of o Cet, and we 
chose to model R Leo with a model identical to the o Cet model, 
but with a luminosity reduced to match the smaller period (310 
versus 332 day s). Miras typically also h ave reduced masses at re- 
duced periods ( I Jura & Kleinmannlll992h . but we chose not to re- 
duce the mass so as to see the differential influence of luminosity 
alone (Wyatt and Cahn 1983 give a mass of 1.04 M© for R Leo). 
The model has a luminosity derived by assuming R Leo was at a 
distance of 1 10 pc, corresponding to a parallax of 9.1 mas, again 
consistent with the Ivan LeeuwenI OOOVt) value of 14.03±2.65 mas 
within 2-CT. 
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Table 5. Parameters of the extended 2 cycles of the R52 series, including the position of the shock fronts. [To be available online only] 



Model 


Phase 


L 


R 


Toff 


SI 


S2 


S3 


S4 






(Lq) 


(Rp) 


(K) 


(Rp) 


(Rp) 


(Rp) 


(Rp) 


360540 


-0.20 


4442 


0.92 


3411 


2.02 


0.93 






360760 


-0.10 


5000 


0.91 


3535 


1.94 


0.92 






361180 


0.00 


8266 


1.06 


3696 


1.79 


1.11 






361720 


0.10 


7425 


1.24 


3338 


1.57 


1.40 






361860 


0.19 


6491 


1.33 


3118 


=► 


1.71 






361900 


0.32 


4700 


1.36 


2844 




2.10 






361920 


0.46 


3729 


1.25 


2793 




2.38 






361940 


0.54 


2795 


1.15 


2714 




2.51 






362020 


0.60 


2901 


1.06 


2851 




2.58 


1.06 




362200 


0.70 


3837 


0.96 


3209 




2.69 


0.97 




362380 


0.80 


4405 


0.92 


3400 




2.79 


0.94 




362560 


0.90 


5083 


0.91 


3545 




2.81 


0.92 




362780 


1.00 


8236 


1.06 


3704 




2.87 


1.10 




363000 


1.10 


7408 


1.24 


3337 




2.93 


1.37 




363180 


1.20 


6374 


1.33 


3097 




2.89 


1.59 




363380 


1.30 


4897 


1.36 


2870 




2.92 


1.77 




363540 


1.41 


4088 


1.31 


2797 




2.81 


1.91 




363600 


1.49 


3063 


1.24 


2674 




=^ 


1.97 




363700 


1.60 


2567 


1.06 


2767 






2.00 


1.05 


363900 


1.70 


3910 


0.96 


3225 






2.00 


0.97 


364120 


1.80 


4456 


0.92 


3409 






1.94 


0.94 


364380 


1.90 


5076 


0.91 


3544 






1.86 


0.93 


364960 


2.00 


8304 


1.07 


3699 






1.76 


1.12 



Table 6. Parameters of the compact 2 cycles of the R5 2 series, including the position of the shock fronts. [To be available on-line only] 



Model 


Phase 


L 


R 


Tcfl 


SI 


S2 


S3 


S4 






(i©) 


(Rp) 


(K) 


(Rp) 


(Rp) 


(Rp) 


(Rp) 


386260 


2.36 


4640 


1.35 


2836 


2.18 








386280 


2.49 


3813 


1.26 


2796 


2.56 








386320 


2.60 


2569 


1.13 


2680 


2.89 


1.10 






386500 


2.70 


3753 


0.99 


3148 


3.09 


1.00 






386660 


2.80 


4198 


0.93 


3332 


3.27 


0.95 






386840 


2.90 


4680 


0.90 


3486 


3.46 


0.92 






387040 


3.00 


8026 


0.99 


3801 




1.02 






387260 


3.10 


7704 


1.18 


3446 




1.27 






387420 


3.19 


6844 


1.30 


3193 




1.49 






387580 


3.30 


5449 


1.36 


2951 




1.70 






387700 


3.40 


4350 


1.34 


2805 




1.85 






387780 


3.49 


3585 


1.27 


2743 




1.93 






387860 


3.61 


2049 


1.24 


2421 




1.98 


1.09 




388040 


3.70 


3560 


0.99 


3111 




1.97 


0.99 




388240 


3.80 


4311 


0.93 


3355 




1.92 


0.95 




388480 


3.90 


4704 


0.90 


3483 




1.83 


0.93 




388900 


4.00 


8166 


0.99 


3814 




1.69 


1.03 




389540 


4.10 


7818 


1.17 


3473 




1.48 


1.30 




389680 


4.20 


6744 


1.31 


3172 




=► 


1.69 




389720 


4.32 


5120 


1.35 


2908 






2.06 




389740 


4.49 


3819 


1.26 


2803 






2.45 




389780 


4.60 


2700 


1.12 


2728 






2.65 


1.10 


389940 


4.69 


3700 


0.99 


3128 






2.77 


1.00 


390120 


4.80 


4182 


0.93 


3332 






2.88 


0.95 


390300 


4.90 


4674 


0.90 


3482 






3.01 


0.93 



It is clear in Figure [6] that this model is too blue to be an ef- 
fective model for R Leo in both V-K and J-K colours, and has 
a visual amplitude much larger and a K amplitude slightly larger 
than R Leo. The amplitude of the R Leo model is slightly smaller 
than the amplitude of the o Cet model due to the reduced luminos- 
ity. This luminosity change by itself is not enough to explain the 
different visual amplitudes of the real R Leo and o Cet. It is pos- 



sible that these two stars have different metallicity (i.e. the redder 
colours of R Leo could be because it has a higher metallicity). A 
difference in mass is also possible, with a compensating change in 
luminosity within that allowed by the parallax error to retain the 
same period. 

Like the o54 series, we computed an extended sub-series (day 
numbers 3600-4500 in Figure |2j and a compact sub-series (day 
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Table 7. Parameters of the C5 series, including the position of the shock fronts. [To be available on-line only] 



Model 


Phase 


L 


R 


Teff 


SI 


S2 


S3 


S4 S5 






(i©) 


(Rp) 


(K) 


(Rp) 


(Rp) 


(Rp) 


(Rp) (Rp) 


375360 


-1.30 


2673 


0.92 


2532 


2.49 


0.92 






375490 


-1.20 


3407 


0.85 


2810 


2.56 


0.87 






375630 


-1.10 


3529 


0.80 


2924 


2.59 


0.83 






375780 


-1.00 


6640 


0.87 


3271 


2.63 


0.92 






375920 


-0.90 


7218 


1.03 


3077 


2.62 


1.13 






376050 


-0.80 


7039 


1.16 


2883 


2.60 


1.37 






376160 


-0.70 


6412 


1.23 


2738 


2.55 


1.57 






376260 


-0.60 


5537 


1.24 


2627 


2.48 


1.74 






376330 


-0.50 


4263 


1.21 


2493 


2.38 


1.84 






376380 


-0.40 


2851 


1.10 


2357 


2.19 


1.91 






376470 


-0.30 


2678 


0.92 


2543 


=> 


1.93 


0.92 




376630 


-0.20 


3392 


0.84 


2816 




1.97 


0.87 




376820 


-0.10 


3595 


0.80 


2937 




1.97 


0.82 




377110 


0.00 


6722 


0.88 


3261 




1.86 


0.94 




377490 


0.10 


7253 


1.04 


3058 




1.75 


1.18 




377710 


0.20 


7037 


1.16 


2876 




=> 


1.53 




377750 


0.30 


6400 


1.23 


2735 






2.01 




377760 


0.40 


5465 


1.24 


2618 






2.39 




377770 


0.53 


3720 


1.21 


2407 






2.83 




377790 


0.61 


2633 


1.18 


2233 






3.06 




377880 


0.70 


2883 


0.91 


2606 






3.31 


0.91 


378020 


0.80 


3370 


0.84 


2823 






3.55 


0.86 


378160 


0.90 


3813 


0.80 


2981 






3.76 


0.82 


378320 


1.00 


6750 


0.89 


3253 






3.96 


0.93 


378470 


1.10 


7246 


1.06 


3035 






4.15 


1.18 


378570 


1.21 


6955 


1.18 


2852 






4.34 


1.44 


378630 


1.31 


6250 


1.23 


2713 






4.52 


1.65 


378660 


1.40 


5395 


1.24 


2610 






(4.66) 


1.80 


378680 


1.49 


4166 


1.21 


2476 






(4.80) 


1.91 


378710 


1.61 


2512 


1.06 


2326 






(4.97) 


1.97 


378810 


1.70 


2972 


0.90 


2631 








1.95 0.91 


378960 


1.80 


3380 


0.83 


2827 








1.89 0.86 



Table 3. Pai'ameters of the compact o5 4 cycle, including the position of the 
shock fronts. [To be available online only] 



Model 


Phase 


L 


R 


T.ff 


SI 


S2 


S3 






(Lq) 


(Rp) 


(K) 


(Rp) 


(Rp) 


(Rp) 


248480 


-0.30 


3243 


0.99 


2982 


2.03 


0.99 




248680 


-0.20 


4122 


0.90 


3312 


2.02 


0.91 




248900 


-0.10 


4379 


0.86 


3450 


1.99 


0.87 




249240 


0.00 


7870 


0.95 


3793 


1.92 


0.97 




249960 


0.10 


8358 


1.18 


3458 


1.80 


1.29 




250360 


0.21 


7432 


1.34 


3152 


^ 


1.69 




250400 


0.29 


6239 


1.39 


2960 




2.10 




250420 


0.38 


4768 


1.41 


2749 




2.42 




250440 


0.53 


3520 


1.33 


2628 




2.96 




250460 


0.60 


2847 


1.32 


2497 




3.18 


1.13 


250640 


0.70 


3305 


0.99 


2990 




3.46 


0.99 


250820 


0.80 


4019 


0.91 


3287 




3.74 


0.91 


250980 


0.90 


4238 


0.85 


3436 




4.01 


0.86 


251160 


1.00 


7713 


0.96 


3768 




4.42 


0.97 



numbers 8200-9100 in Figure|2ll. The compact sub-series had its 
detailed radiative transfer model truncated at 4 i?p, because the very 
low density in the outer layers (< 10^^^ gcm^'^) had extremely 
low opacity and the models had outer-layer physical conditions out- 
side the range where our equation of state was valid. Again, the V, 
J and K fluxes were similar in each case (only the extended sub- 
series is displayed). 



2.3 C81 8160 L© series 

With the same input physics as the o Cet models, we attempted to 
create a longer period series appropriate for the Mira variable R Cas 
(P = 430 days) by increasing the luminosity and mass. However, 
the model pulsation amplitude became much too large, requiring us 
to increase the ai, parame ter. The HIPPARCOS distance for R Cas 
in lWhitelock et alj (2003) would require the star to be very under- 
luminous when compared w ith the mean solar-vicinity P-L rela- 
tion ship of Whitelock et al] ( |2008|) or the LMC P-L relationship 
(e.g Hughes & Woodlll990D . For this model series, the luminosity 
of our model was derived by as suming R Cas falls on the mean 
solar-vicinity P-L relationship of IWhitelock et al.l J2008h and is at 
a distance of 204 p c, with a corre s pondi ng parallax of 4.9 mas, al- 
most 3-(T from the Ivan LeeuwenI ( 120071) value of 7.95±1.02mas. 



The mass adopted for R Cas is 1.35 M© as given bv lWvatt & CahnI 
[1983) . 

The light curves of R Cas in the V, J and K bands are shown 
in FigurejS] corresponding to day numbers ~4000-5700 chosen for 
detailed radiative transfer computation in Figure [3] The J and K 
light curves of the model fit the observations quite well, while in all 
cycles where the near-maximum model was computed, the models 
are too bright in V near maximum light. The near-maximum con- 
tinuum effective temperature is ~3800K in Table [8] corresponds 
to an MO or Ml g iant according to the temperature calibration of 
iFluks et alj ( 119941) and appears like an M2 giant in the TiO fea- 
tures as predicted by our model spectra. This is much too warm for 
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Table 8. Parameters of the C81 series, including the position of the shock fronts. [To be available on-line only] 



Model 


Phase 


L 


R 


Teff- SI 


S2 


S3 


S4 


S5 S6 






(i©) 


iRp) 


(K) (fip) 


(i?p) 


iRp) 


iRp) 


{Rp) (Rp) 


243570 


-0.40 


3917 


1.58 


2182 2.74 


1.96 








243670 


-0.30 


3705 


1.34 


2336 2.40 


2.06 


0.98 






244000 


-0.20 


5772 


0.87 


3235 =^ 


2.14 


0.87 






244140 


-0.10 


5133 


0.79 


3311 


2.29 


0.79 






244330 


0.00 


9462 


0.81 


3811 


2.25 


0.81 






244700 


0.10 


12460 


1.08 


3522 


2.21 


1.14 






246230 


0.20 


12120 


1.28 


3212 


2.20 


1.49 






246810 


0.30 


11150 


1.40 


3015 


2.05 


1.87 






246910 


0.40 


8327 


1.44 


2758 


=> 


2.31 






246940 


0.52 


5324 


1.65 


2310 




2.82 






246960 


0.61 


4427 


1.84 


2087 




3.07 






247020 


0.70 


4103 


1.63 


2175 




3.39 


1.01 




247280 


0.80 


5302 


0.89 


3130 




3.67 


0.89 




247450 


0.90 


5184 


0.79 


3322 




3.88 


0.79 




247560 


1.00 


10470 








4.13 


0.78 




247690 


1.10 


11640 


1.06 


3498 




4.39 


1.08 




247810 


1.20 


12060 


1.27 


3228 




4.54 


1.35 




247900 


1.30 


10990 


1.39 


3014 




(4.80) 


1.60 




247970 


1.40 


8416 


1.44 


2769 




(4.92) 


1.79 




248020 


1.50 


5748 


1.58 


2405 






1.91 




248060 


1.59 


4234 


1.67 


2169 






1.99 




248130 


1.70 


4143 


1.41 


2341 






2.00 


1.00 


248430 


1.80 


5493 


0.88 


3175 






1.96 


0.88 


248580 


1.90 


5250 


0.79 


3326 






1.89 


0.79 


248730 


2.00 


9757 


0.79 


3882 






1.76 


0.79 


249030 


2.10 


12020 


1.06 


3523 






1.59 


1.09 


249400 


2.20 


12350 


1.26 


3255 






=> 


1.41 


249500 


2.30 


11070 


1.38 


3029 








1.82 


249580 


2.39 


8557 


1.42 


2801 








2.09 


249710 


2.50 


5655 


1.57 


2400 








2.35 


249890 


2.60 


4468 


1.68 


2191 








2.53 


250020 


2.70 


3955 


1.59 


2180 








2.68 


250290 


2.80 


5609 


0.89 


3178 








2.79 


250440 


2.89 


5270 


0.80 


3311 








2.88 


250550 


3.00 


9899 












2.94 


250670 


3.10 


11740 


1.06 


3501 








2.97 


250780 


3.20 


11920 


1.27 


3217 








2.98 


250870 


3.30 


10810 


1.38 


3012 








2.98 


250960 


3.40 


8299 


1.43 


2765 








2.95 


251040 


3.49 


5758 


1.55 


2429 








2.87 



R Cas which has a catalogued spectral type of M6-M10. Therefore, 
we are forced to conclude that Teff is too high for this model. 

2.4 C 5 5050 Lq series 

As the 8160 L© model or R Cas was so clearly discrepant near max- 
imum, and as the individual HIPPARCOS distance in Whitelock et 
al. (2008) would give a luminosity of only ~3770 L© , we chose 
to construct a lower luminosity model. The luminosity was deter- 
mined by fixing the mixing length parameter which was decreased 
to a more standard value of Om ~ 2. The luminosity was then tuned 
to match the model period to that of R Cas. As usual, the turbulent 
viscosity parameter was then tuned to match the bolometric ampli- 
tude of R Cas and the model. The resulting luminosity (5050 L©) 
suggests that R Cas is at a distance of 166 pc with a corresponding 
parallax of 6.0 mas: this is now consistent with the HIPPARCOS 
value within 2-a. 

The near-maximum effective temperature of this model is now 
~ 3250 K, both consistent w ith an M6 spectral type from the cal- 
ibration of iFluks et alj J 19941) and providing a better match to the 



observed M6 spectra in iFluks et"ai] ( 1 19941) than either M5 of M7 
spectra. It also has a a V-K colour that matches that of R Cas (see 
Figure |71l. 



2.5 Using the Models 

For each model phase in Tables [3] through [S] we provide the 
full model output q Indeed, this output was already used by 
IWoodruff et al.l ( 120090 in advance of publication in order to com- 
pare the models to wavelength-dispersed infrared interferometry. A 
sample 3 lines from one of these tables is given in Table|9] In order 
to use these models to compare to a specific observation, integra- 
tion over a filter profile F(A) is required: 



Lf ^ L{\)dX 



http://www.physics.mq.edu.au/~mireland/codex/ 



(2) 
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Table 9. An excerpt of the table available on-line describing the model output for the o54 series model 285180. The center-to-limb variation (CLV) is shown 
as a fraction of the central intensity /q as a function of normalized radius on the apparent stellar disk, r/5Rp . 



Wavelength L^ 



lo 



CLV: r/5Ri, 



and 



I/Io 



(^m) 



(ergs^fm-i) (erg s cm'^ /^m"! sr"!) 0.020 0.100 0.140 0.170 0.190 0.220 0.290 0.360 0.440 



1.598 0.9786E+37 0.1124E-I-10 

1.599 0.8160E+37 0.3054E-I-09 

1.600 0.1106E+38 0.1427E-I-10 



0.997 0.923 0.828 0.691 0.451 0.082 0.066 0.053 0.046 
0.999 0.977 0.955 0.932 0.913 0.879 0.796 0.662 0.337 
0.998 0.940 0.863 0.743 0.274 0.024 0.022 0.022 0.026 



M=1.35Mq L=5050Lg Z = O.Oa P = 427 days 




2000 



4000 6000 

time(days) 



8000 



10* 



Figure 4. The luminosity and mass-zone positions in the C50 series, show- 
ing the locations of the phases chosen for detailed model computation. 
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Figure 5. Light curves of the o 5 4 series m J, K and V bands compared 
with observations of o Cet bv .Whitelock et al. C2000,) and the AFOEV. 
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Figure 6. Light curves of the R5 2 series i n the extended cycle in J , K and 
V bands compared with observations by Iwhitelock et alj 1200(1) and the 
AFOEV. 





Figure 7. Light curves of the C5 seri es J, K and V bands compared with 
observations bv lNadzhip et alj J200lh and the AFOEV. 
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Figure 8. Light curves of the C8 1 seri es J, K and V bands compared with 
observations bv lNadzhip et alJJ200il) and the AFOEV. 
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Figure 9. A H-band comparison of the R Cha spectrum at phases 0. 1 (top) 
and 0.3 (bottom) to the 287820 (top) and 250460 (bottom) models of the 
o54 series. 



Figure 10. A K-band comparison of the R Cha spectrum at phases 0. 1 (top) 
and 0.3 (bottom) to the 287820 (top) and 250460 (bottom) models of the 
o5 4 series. 



Ipix) 



f F{X)Io{\)fix,\)d\. 



(3) 



Here Lf is the stellar luminosity as seen through the filter, 
If{x) is the intensity profile seen through the filter, Io(A) are the 
tabulated values of the central intensity and f{x, A) are the tabu- 
lated values of the normalized CLVs. It is much more preferable 
to use realistic filter profiles F[\) with smooth edges (e.g. a Gaus- 
sian) rather than square-edged filters in order to minimise noise due 
to the opacity sampling. Interferometric visibilities can then be ob- 
tained from the Hankel transform of If{x). 



3 MODEL SPECTRA 

The spectra computed in the CODEX models using the default 
wavelength table come from an opacity sampling method with a 
spectral resolution of up to ~ 10^. However, in order to accurately 
compare with observations at any wavelength, at least ~100 wave- 
lengths have to be averaged together, preferably using a non-square 
edged filter. In turn, this means that these default model outputs can 
only be used at a spectral resolution of R ^ 100 or lower. This is 
especially true where the CO bands in H and K bands are con- 
cerned, where there is a combination of very strong absorption and 
near-continuum spectral features. 

For the purposes of comparing spectra with observations in 
H- and K-bands, we therefore also computed spectra and center-to- 
limb variations at R~ 10^ in these bands. In Figures [9l and 1 101 we 
convolved the model spectra with a Gaussian of Full-Width Half 
Maximum equivalent to a spectral resolution of R — 1000, and 
compared t he model spectra with ob servations of the Mira variable 
R Cha from lLancon & WoodfeOOOl) . We chose to compare with the 
o54 series, because the parameters of R Cha are most like those of 
o Cet. There was an arbitrary scaling applied to the observed spec- 
tra. These factors differed by 0.05 mags between the H- and K- 
bands for the phase 0. 1 spectrum and 0. 15 magnitudes for the phase 
0.3 spectrum: these small differences . The best model fits were 
from phases 0.3 and 0.6, where temperatures were ~500 K cooler 
than phases applicable to the R Cha observations. This demon- 
strates that in the CO overtone bands, the models are too hot for 
R Cha. To make this statement more quantitatively based on spec- 
tral synthesis, metallicity effects would have to be considered also 
(beyond the scope of this paper). 



4 OBSERVATION AND MODEL COMPARISON: O CET 

The grand total of all available observations for Mira variables with 
parameters similar to those of the model series presented here is far 
too vast to compare to the model series of this paper in a concise 
manner. Therefore, we have chosen to examine the available ob- 
servations of o Cet in a general sense to describe the similarities 
and differences between the o54 series and o Cet to further illus- 
trate the utility of the model series and the wealth of information 
available to constrain models. 

Table \T0\ summarizes most of the key observations available 
for o Cet. Time-dependent photometry is available between ultra- 
violet and radio wavelengths, with the best light curves available in 
the V, J, H and K bands, as shown in Figure [5] (V,J,K). There is 
reasonable agreement between the model and observations for light 
curve shape, amplitude and visible-infrared phase offsets. 

Spectral classification should also give an observed effec- 
tive temperature. U nfortunately, the MK spectral classification (e.g 
iKeenan et al J 1 19740 is based on B and V bands, w here non-LTE 
effect s in an extended atmosphere are very strong Jlreland et al.l 
I2OO81) . Consequently, the combination of effective temperature and 
metallicity can not be directly fit to observations. Spectra of o Cet 
are also not available electronically to the knowledge of the au- 
thors - a modern library of bright Mira spectra would certainly be 
of great use to future modelling efforts. In particular, infrared spec- 
tra are a much more reliable model output, and phase-dependent 
infrared spectra would be a wonderful tool for tuning model pa- 
rameters. 

Resolved observations have been made at wavelengths be- 
tween 346 nm and 7 mm, with the broad range of highly wavelength 
dependent diameters being consistent with Models. The for m of an- 
gular diameter versus wavelength curve as shown in lWoodruff et al.l 
(.2009.) between 1 and 4 microns was very similar to the o5 4 model 
series, however, in that paper, the models were placed at a distance 
that best fit the angular diameters. If instead the models are placed 
at the distance that best fits the K-band photometry, the angular di- 
ameters as a function of phase are given in Figures [TT] to [T3] It is 
clear that the mean diameter of the models are too small, and that 
the phase-dependence of the observed diameter is less pronounced 
than in the model. Possible solutions to this are given in Sections|5] 
and|7] 

In the radio, the measured angular size of o Cet corresponds 
to the angular size of the Na and/or K ionisation edges. The con- 
sistency between models and observations shows that local thermal 
equilibrium is a reasonably assumption for defining the ionization 
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Figure 11. Model diameters of the o54 series in a narrow 1.24 ^m band- 
pass based on fitting to a single spatial frequenc y where V = 0.6, with the 
measured diameters from lWoodruff et alj J2008l) over-plotted. 




Figure 12. The same as Figure fTTI except for the H-band. 



fraction. In the ultraviolet, the measured FWHM of 35 mas, corre- 
sponding to a shell uniform disk diameter of ~56 mas, will pro- 
vide a strong constraint on the radii of small dust grains. However, 
model outputs are not currently available for that wavelength range. 



5 PREDICTING FUNDAMENTAL PARAMETERS WITH 
MODELS AND OBSERVATIONS 

For each Mira variable modeled in this paper, there are currently 
3 physical (M, L and Z) parameters and 2 model parameters {a,n 
and a^). In the parameter neighborhood of the o54 series, we find 
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Figure 13. The same as Figure [171 except for the K-band. 



that the radius of the "parent" star approximated by linear pulsation 
is given by: 



R* 



L, 



M* 



2167?. 



'0 



^ 5400L, 



© 



-1.1M( 



or in angular units: 



oc d'-'F^^-^M-^-^a, 







-0.7/ Z .0,2 



^3.5^ ^0.02^ 



(4) 



(5) 



where d is the distance, and F is the received wavelength- 
integrated stellar flux. In this section we will only discuss the most 
direct measurements of effective temperature, derived from near- 
continuum interferometry and photometry, as spe ctral fitting in the 
presence of non-LTE effects ilreland et al.ll2008i ). metallicity and 
abundance errors has yet to be demonstrated for extended M giants. 

There is no dependence of the radius on a^/ , but there is a small 
dependence of the period P on Qi, which we will neglect here. The 
period, which typically is expressed in terms of mass and radius 
(the so-called PMR relationship), we will express in terms of our 
model parameters L,, M, and am' 



L, 



M. 



261d 



' 5400L 







I.IM, 







-l.S/ Qm 

4.5 



)" (— ^) 



(6) 



Although these relationships are only approximate and do not 
hold over a wide range of parameters, they demonstrate the com- 
plex interplay between the model input parameters. In principle, a 
measurement of Z from spectral synthesis, and measurements of 
period, amplitude, angular diameter, luminosity and distance are 
enough to constrain M, am and a^. However, a 10% distance un- 
certainty (the best of any nearby Mira) translates into a 20% L un- 
certainty, or a 20% mass uncertainty keeping everything else fixed 
at a given period. The relationships are further complicated by the 
non-linear pulsation period differing significantly from the linear 
pulsation period (e.g. Table 1), depending on amplitude. 

Consider first the problem of calibrating mass independently 
of pulsation models. Orbital periods for non-interacting Miras are 
far too long for combined visual and spectroscopic combined orbits 
to obtain dynamical masses due to the large radii of Miras. Clusters 
form a potential hunting ground for Mira variables where the AGB 
can b e calibrated at a known initial mass (e.g. Lebzelter & Wooq 
I2007D . but the final mass of the Mira is a function of the assumed 
mass loss history, clusters do not easily provide the same age and 
metallicity range of Miras as in the field, and a direct radius mea- 
surement is not yet possible. Mass can also be estimated from kine- 
matics. This is best done for Miras with kinematics inconsistent 
with the thick disk or halo. The best example of this for nearby 
Miras is o Cet. 

o Cet has a {Uy,W) space velocity of (-26, -6 2, -89) km s~^ 
when using the revised HIPPARCOS distance from van LeeuwenI 
( 120070 . This space velocity is unusually large for a Mira, which is a 
major reason why the interaction between o Cet and the inters tellar 
medium produces such an impressive tail jMartin et alJuOOT ). Al- 
though kinematics is often inconclusive when applied to individual 
stars, this space velocity falls within t he 97% probability contour 
for thick disk memb ership accor d ing tolReddv et alj J2006l) . In the 
detailed analysis of iRobinet alj yOOSf), the W velocity of o Cet 
is inconsistent at 5g with even the old (5-10 Gyr) thin disk, and is 
most consistent with being a member of the thick disk, modeled as 
a single stellar populat ion of 11 Gyr age. Imp ortantly, o Cet can not 
be a runaway star (e.g. lHoogerwerf et al.l200 0) where its space ve- 
locity is due to a single strong gravitational interaction in its past, 
because it has retained its wide companion Mira B, and the orbit 
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Table 10. Summary of key observational data on o Cet. Units: Vega magnitudes for light-curves, mjy for radio photometry, milli-arcsec for diameters. 



Wavelength 


Data Type 


Observation 


Model range 


Phase 


References 


(/im) 




Range 




Coverage? 




0.55 


Light Curve 


2.5-9.5 


2.0-1 1 


Y 


AAVSO,AFOEV 


1.2 


Light Curve 
Light Curve 


(-1.7)-(-0.3) 
(-2.6)-(-1.3) 


(-1.8)-(-0.6) 
(-2.6)-(-1.5) 


Y 
Y 


Whitelock et al. QOOO) 


1.65 


Whitelock et al. (2000) 


2.2 


Light Curve 


(-3.0)-(-1.8) 


(-3.0)-(-2.0)'= 


Y 


Whitelock et al. (2000) 


0.4-0.6 


Spectral Type 


M5e-M9e 


None'' 


Y 


Skiff (2009); Samus et al. (2004) 


0.307 


Photometry 


14.85 


None"^ 


N 


jCarovskaetal. (1997) 


0.346 


Diameter 


35 (FWHM) 


None*" 


N 


Karovskaetal. (1997) 


0.45-1.03 


Diameter 


31-103 


25-60^* 


N 


Labeyrie et al. (197T) 


0.68-0.92 


Diameter 


20-60 FWHM 


16-38 


N 


Ireland et al. (2004) 


1.24 


Diameter 


22-30 


17-34 


Y 


Woodruff et al. (2008J 


2.26 


Diameter 


31-37 


17-29 


Y 


Woodruff et al. (2008) 


1.24 


Diameter 
Diameter 
Diameter 


22-30 
25-68 
46-55 


17-32 
17-58 
21-72 


Y 

N 
Y 


Woodruff etal. (2008) 


1.1-3.8 


Woodruff et al. (2009) 


11.15 


Weiner et al. (2003) 


3.6 cm'" 


Photometry 


0.20-0.37 mJy 


0.09-0.32 mJy 


N 


Reid & Menten (1997): Matthews & Karovska (2006) 


7 mm 


Diameter 


52 


29-55 


N 


Reid & Menten (2007) 







" Shorter- wavelength radio observations are consistent with a u^ power law within eiTors. A 0.12 mJy contribution from Mira B has been subtracted. 

''No model spectral type calculations are possible. See text. 

'^Distance to model fixed to 107 pc so that K-band maximum agrees. 

''When the fluorescence scattering approximation is used (as in Paper I), the range becomes 26-80 mas, with the upper diameter limited by the model 5 Rp 

surface. 
■^Although the ultraviolet is important for deep-atmosphere temperature profiles, no ultraviolet fluxes are output because the near-surface opaciites are likely 

unreliable. 



can not be highly eccentric as the periastron must be outside the at- 
mosphere of o Cet. Therefore the pr ogenitor mass of o Cet is almost 
certainly less than 1.1 Mq (e.g. the lOirardi et alj ( I200Q) evolution- 
ary tracks places a 1.06 M© star on the TP-AGB after lOGyr, or 
less for sub-solar metallicity). This places the current mass of o Cet 
at 1 Mq or less. This is close to our model value of 1.1 A4q. How- 
ever, if L were to be kept constant and A/ decreased, am would 
have to be increased further from its already large value in order to 
maintain the period. We will discuss this further below after first 
discussing Um- 

The comparisons of models to observations in this paper has 
already provided significant evidence that model temperatures are 
too high, implying that Um = 3.5 is too high a value for Miras with 
parameter ranges applicable to o Cet and R Cas. There are 3 key 
pieces of evidence: the near-maximum V-K colours, the infrared 
spectral fitting for R Cha, and the measured angular diameters. 

All series with Um = 3.5 are too warm near maximum, 
as shown by their V-K colours of <5.0. The R52 and C81 se- 
ries compared to R Leo and R Cas respectively are the most dis- 
crepant. Although in general this could be due to fundamental pa- 
rameters M and L only, in this case it is not possible. The RCas 
R81 1.35 Mq, 8160 Lq model is reasonable for a Mira: i.e. all 
AGB stars above a certain mass will go through an 8160 Lq phase, 
and the mass of 1.35 M© at 8160 Lq will result from some ini- 
tial mass/rhe_systematic study of Mira spectral type at maximum 
bv lKeenan et al.l ( 119741) found a trend of systematically later spec- 
tral types with increasing Mira period, with no Miras having near- 
maximum spectral types earlier than M5 in the 350-500 day pe- 
riod range. The model spectral type of ~M2 for R Cas is therefore 
much too warm. This means that the combined choice of Um = 3.5 
and ai, = 0.32 together is incorrect for stars in the vicinity of 
M = 1.35 Af© andL = 8160 Lq. 

R Cha is a Mira variable very much like o Cet (period, ampli- 
tude, colour), except for its more uncertain distance. In Section[3] 



we saw that the models were ~500 K too hot for R Cha. However, 
R Cha has the hottest near-maximum spectral type of all M ira vari- 
ables at periods of ~300-350 days in lKeenan et al.l ( 1 1974) . so this 
suggests that the models are too hot for any Mira. Inspection of 
some low-metallicity test models (to be published) indicates that 
decreasing the metallicity reduces this discrepancy but does not re- 
move it. 

Although the range of model diameters reported by 
I Woodruff et al.l (|2008J were roughly consistent with the range of 
diameters for the Miras studied in that paper, the minimum mea- 
sured diameter for each star was larger than the minimum diam- 
eter predicted by models. This remains true for the current model 
series (see Section|4ll, especially for o Cet. A discrepancy in near- 
continuum diameters at phases ~0.9-0.2, where the photosphere is 
relatively compact can only be rectified with an increase in model 
radius by 10-20%. 

Let us consider how models of o Cet could be modified in 
order to achieve a larger apparent radius. The simplest method is 
to place the star at a closer distance with a lower luminosity and 
the same linear radius and mass (preserving the period). We note 
that the most recent HIPPARCOS analysis of Ivan Leeuwerj J2007l) 
places o Cet as close as 75 pc within 2a, which would mean lumi- 
nosities as low as 2650 Lq are consistent with direct observations. 
Such a low luminosity would, however, be 0.8 magnitudes below 
the LMC P-L relationship (Whitelock et al. 2008). This relation- 
ship provides a stronger constraint on the absolute K-magnitude of 
o Cet, as the dispe rsion in the relationship is only 0.13 magnitudes 
dFeast et al.lll989l) . Adopting this relationship places a 2a lower 
limit on the o Cet luminosity of 4400 Lq . By solving the approxi- 
mate equations|4]through[6]after decreasing the Al to 1 Mq and in- 
creasing angular diameter by 10%, the model luminosity becomes 
4200 Lq and the mixing length parameter am = 2.9 (effects of ai, 
and Z are neglected here). It is therefore likely not possible to fit all 
observations of o Cet unless the model mass is greater than 1 Mq 
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Figure 14. Temperature (K ; upper left), density (g cm ^ ; lower left panel) and the moecular density Pmol for molecules H2O (upper right) and TiO (lower 
right) as a function of r/Rp for 3 phases of a cycle of the o54 model series: 285180 (phase -0.20, solid Une), 285860 (+0.20, dashed), 286100 (+0.49, 
short-dashed). 



or the luminosity more than 2a below the LMC P-L relationship 
- each of which would be a controversial claim. This discussion 
shows just how difficult it is to provide precisely calibrated models 
of Mira variables without clear observational reference points. 

The key problem here of finding the best value for am as a 
function of M and L can be expressed as a problem of finding the 
radii of real stars at given M and L. Mixing length only provides a 
way to calculate R for given M and L once free parameters (espe- 
cially the mixing length) are fixed based on known stars. Models of 
main-sequence stars and even K-giants can be calibrated very well 
from the sun and other stars of accurately known parallax. How- 
ever, extrapolating to M giants, where pressure scale-heights can 
be a significant fraction of a radius, is not expected to be reliable. 

What is therefore needed is a sample of M giants of well- 
known M, L and R, from which to calibrate the mixing length. Al- 
though we can not measure M, period can serve as a proxy for M 
given L and R. L and R can be measured accurately for a sample of 
stars with relatively compact atmospheres (i.e. with a well-defined 
R) that have a well known parallax, photometry and effective tem- 
perature. At this point, spectral synthesis is not reliable enough or 
calibrated well enough for M giants to produce accurate effective 
temperatures and compositions. Therefore, this sample should in- 
clude semi-regular pulsators with accurate photometry, angular di- 
ameters and periods. Either accurate periods or angular diameters 
are currently missing for many of the closest M giants, so we sug- 
gest that measuring and collating such information should be an 
active area of research. 



Assuming that a reasonable value for am can be prescribed 
for Mira models, the main free parameters for any individual Mira 
are composition and a^. Composition (primordial metallicity and 
C/O ratio) must come from spectral synthesis. Given a set of Miras 
with well-measured distances, L is determined from observations 
and M from the period of the Mira - only q„ can be used to tune 
the model amplitude and should be relatively easy to calibrate. 



6 THE MOLECULAR-SHELL SCENARIO 

Interferometric measurements of Mira variables have been shown 
to be internally consistent only if there is a layer of molecular water 
far above the continuum-forming photosphere (e.g. Weiner 20041). 
Observational comparison with models that include this kind of wa- 
ter and/or dust layer have so-far been dominated by non-physical 
models, i.e. those that neither provide a mechanism for elevating 
the emitting material nor calculations of the chemistry that deter- 
mines which components dominate the radiative transfer at which 
radii. 

Nevertheless, these ad-hoc models have provided a relatively 
simple picture for the regions around Mira variables and have im- 
pressively fitted a limited selection of observable properties. In this 
section we will examine the physical and observable properties of 
the molecular shells in our model series. 

Figs. 1 to 4 show that cycle-to-cycle variations in the pulsation 
models are generally quite modest in terms of the luminosity, of the 
temperature at and the position of the Tg=2/3 layer and of the posi- 
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Figure 15. The same as Figure [T4l for the same phase 0.7 of 3 models of different cycles of the o 5 4 model series: 248480 ("extended" series -0.30, solid line), 
250640 ("extended" series +0.70, dashed), 288140 (+1.70, short-dashed). 



tion of deep-layer mass zones below the Tg=2/3 layer. The o54 and 
C81 model series show somewhat more pronounced cycle-to-cycle 
effects than the r52 and C50 models. Inspection of specific numer- 
ical values of the radius R of the rRoss=l layer and the therefrom 
derived effective temperature Tcb oc (L/_R^)^'* given in Tabs. 2 
to 8 for the non-grey atmospheric stratifications confirm this cycle 
stability. 

In contrast, we notice substantial differences between different 
cycles, and often between successive cycles, in terms of the posi- 
tions of high-layer mass zones (Figs. 1 to 4). These differences are 
closely related to substantial differences of the strengths and posi- 
tions of outward traveling shock fronts. Inspection of shock-front 
positions (Tabs. 2 to 8) in the selected cycles for which detailed 
atmospheric models were computed show a shock front typically 
emerging at pre-maximum phase around -0.3 to -0. 1 , then traveling 
outward during about 1 to 1 1/2 cycles while it becomes weaker 
and slower before the subsequent front catches up and both fronts 
merge. Typically, the outer front starts retreating before merger, 
but occasional shock fronts traveling towards circumstellar space 
(where they eventually fade away) are also seen in the Tables. 

The positions and heights of shock fronts at different phases 
and in different cycles determine the upper atmospheric density 
stratification and, therefrom, the details of the temperature strati- 
fication and of the partial-pressure stratification of molecular ab- 
sorbers. The assumptions of local thermodynamic equilibrium and 
spherical symmetry are sufficient to derive this stratication (as may 
not be the case for the so-called MOLsphere ( Tsuji 2000) in super- 
giants, due to the co-existence of the chromosphere in those stars). 



The st udy of IXei et alj ( 120030 . based on models of iHofmann et al.l 
( Il998[) . shows that the details of shock-front propagation may lead 
to strong cycle-to-cycle differences of the stratification of the outer 
atmosphere resulting in strong differences of the density and geo- 
metric characteristics of water "shells", i.e. of layers whose absorp- 
tion is dominated by water molecules. 

Figures[T4]to[T6]demonstrate, for the o54 and the r50 series 
of the here presented CODEX model sets, the drastic phase and cy- 
cle effects of shock-front propagation on the temperature-density 
stratification and on the appearance of H2O and TiO "shells". Den- 
sity decreases monotonically with radius, while the sharp decrease 
at a shock front provides the outer edge of any "shell". The de- 
crease of temperature with radius provides a relatively sharp edge 
to the region where water can exist in chemical equilibrium, and 
this provides the inner edge of any "shell". At most infrared wave- 
lengths, gas is reasonably transparent between temperatures where 
H- opacity is dominant (> 3000 K) and where water is domi- 
nant (< 1800 K). A similar pattern is found for the TiO molecule 
which, however, is formed in somewhat deeper layers than water 
and, therefore, does not depend so strongly on upper-atmosphere 
shock-fronts and shows smaller, though by no means negligible, 
cycle-to-cycle effects. 

The model-predicted effects of such molecular "shells" upon 
the absorption properties of the stellar atmosphere at observation- 
ally important wavelengths have bee n discussed and compa red to 
observations by llreland et all ( l2008h , IWoodruff et all ( [2009,) and 
Wittkowski et al. (in preparation). Figure [TT] shows typical cycle- 
to-cycle differences seen in the water-contaminated spectrum of the 
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Figure 16. The same as Figure [74l for the same phase 0.0 of 3 models of successive cycles of the r50 model series: 375780 (-1.00, solid line), 377110 (0.00, 
dashed), 3783200 (+1.00, short-dashed). 



3 models presented in Figure[T5] Typical effects of water "shells" 
upon the shape of the centre-to-limb variation have been discussed 
in the model study of lTei et alj J2003h . 

Though semi-empirical "shell" scenarios have been used with 
remarkable success for interpreting spectroscopic a nd interferomet- 
ric o bservations of absorption features of H2O (e.g Matsuura e t al.l 
200a: Mennesson et al. 2002; Ohnaka 2004; Perrin et al. 20041; 



Weineil l2004h and of TiO (e.g iReid & Goldstorj l2002h . models 



show that such scenarios can at best provide a very rough picture 
of the approximate instantaneous position and extent of absorbing 
layers. Such semi-empirical models cannot provide any informa- 
tion on changes of these layers with phase (since local molecule 
abundance and resulting molecule absorption as a function of local 
values of p(r) and T{r) sensitively depends on details of shock- 
front progression). Note also that absorption by CO in low exci- 
tation lines extends from the continuum-forming photosphere at 
^-^3000 K right to the w ind region, so a CO "shell" scenario (e.g. 
iMennesson et al.ll2002r) should be considered with particular cau- 
tion. 
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Figure 17. Spectra of the 3 same-phase models of Figure fTH showing cycle- 
to-cycle differences of the spectrum in the J-band region of the spectrum. 



Finally, we note that the existence of shell-like structures no- 
ticeably change the computed spectra in the models, but the effects 
are not so strong that the detailed shell structure can be directly 
inferred from low spectral resolution observations. Modeling high 
spectral resolution observations is beyond the scope of this paper, 
but such a study would have to take i nto account the veloci ty struc- 
ture of the atmosphere explicitly (e.g lNowotnv et alj2010l) . 



7 RADIATIVE ACCELERATION AND MASS LOSS 

As described in lHofnen ( 1200 a) , the conditions for radiative acceler- 
ation to drive mass loss are that dust must be able to form, and that 
the opacity exceeds the critical opacity: 



K> = 12830 cm^g i(_-)(_) \ 

L, Mq Lq 



(7) 
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Figure 18. Interferometric diameters for the o54 series model 286060 
(phase 0.41), fit to the spatial frequency where visibility V = 0.5. The solid 
line corresponds to the default number of dust nuclei, and the solid red line 
corresponds to a facto r of 3.1 incre ase in grain radius, at log(A'^nuc)=-13.7. 
Observations from Wo odruff et al J (2008 ) at phase 0.3 are over-plotted. De- 
spite still not being able to drive a wind, the large diameters at short wave- 
lengths are clearly inconsistent with observations. 



This opacity is attainable at solar metallicity with fully- 
condensed iron-rich du st, which is no t stable until approximately 5 
continuum stellar radii jWoitkduOOm . This opacity is also reached 
by forsterite (i.e. Fe-poor silicate) grains of 400 nm radius due to 
strong scattering (Hofner 2008). 

The CODEX models have a chemical equilibrium model of 
dust formation, so are not appropriate for modeling the slowly- 
growing Fe-rich dust at <5Rp, and indeed we artificially cut-off 
Si condensation at a condensation fraction of 0.25 for this reason. 
However, as discussed in Ireland & Scholz ( 2006), the prescription 
we use for dust formation is reasonably accurate for Fe-poor sili- 
cates. 

The strongest observational constraints on the radii of domi- 
nant dust species are observations that probe the opacity at short 
wavelengths where dust scattering is dominant, and wavelengths 
where water absorption is dominant. As optically thin scattering 
does not affect the spectrum, the best observations to probe this 
difference are resolved observations of Miras as a function of wave- 
length. Figure [TS] shows the interferometri c diameters of o Cet as 
a function of wavelength as measured by I Woodruff et alj J2008I) 
along with the diameters predicted by the 286060 model and the di- 
ameters predicted by the same model with the base- 10 logarithm of 
the number of dust nucleii per H atom TVnuc decreased from -12.2 
to -13.7, and the corresponding maximum grain radius increased 
from 63 to 194 mn. These large grains still have insufficient opacity 
to overcome gravity in this model. We could not increase the grain 
radii further without the optical depths becoming too large at our 
chosen 5 Rp surface, and the preferred value from Hofner ( 2 0081) 
of log(-/Vnuc=-15) would produce diameters that are far too large. 
It is clear that these large grain radii are not consistent with the 
relatively large H2O column densities at several continuum radii 
inferred from infrared interferometric observations. Either radia- 
tion pressure on small Fe-rich grains or large Fe-poor grains could 
drive winds from M-type Mira variables, but the base of the wind 
and the grains that drive it must originate from layers higher than 
where H2O is seen in the near-infrared, i.e. higher than about 3-4 
continuum-forming radii. 



8 CONCLUSIONS AND FUTURE WORK 

The atmospheric models presented here, based on self-excited pul- 
sation models and opacity-sampling treatment of radiation trans- 
port, provide a fairly realistic approximation of the atmospheric 
density-temperature stratification. Many spectral features are pre- 
dicted with satisfactory accuracy but some, like TiO bands, require 
further improvement of the models e.g. the non-LTE treatment of 
Paper I. The present sequence of models also comprises only 4 
combinations of basic stellar parameters at only a single compo- 
sition and, therefore, can only describe a relatively small subset of 
Mira variables. Predictions of the 4 model series presented here 
(Table 1) are available on-line (Section|23](- 

There are several potential causes of substantial deviation 
from spherical asymmetry in Mira variables - inclu ding convectiv e 
cells, weak chaos and Rayleigh-Taylor instabilities IWoitkd ( 120061) . 
Asymmetries are clearly not considered in our model series as 
they are spherically symmetric, but asymmetries are relatively com- 
mon in Mi ra variables when ob served at sufficient angular reso- 
lution (e.g. iRagland et alj 12006)). The one prediction that can be 
made from the models is that the wavelengths most susceptible to 
cycle-to-cycle variations (e.g. L-band where the water shells are 
optically-thick) should also show asymmetries, as high layers on 
opposite sides of the star should not be strongly causally connected 
and show weak chaos. We finally suggest that an observation so far 
missing in the literature is the astrometric motions of the radio pho- 
tosphere over several cycles, which should be a strong indicator for 
the degree of high-layer chaos in Mira atmospheres. 

Determination of the internal fundamental model parameters, 
i.e. mixing-length Um and turbulent viscosity a^, (Table 1), would 
require observation of a set of stars with different mass, luminosity 
and pulsation period. For a given pulsation period, a higher-mass 
star must have a higher luminosity (radius), but a similar effective 
temperature. We suggest that further studies of low-amplitude pul- 
sators with Mira-like periods such as R Dor and W Hya may pro- 
vide the key to tuning the mixing length parameter am of Mira 
model series. These are likely higher-mass stars, but with a similar 
effective temperatures to o Cet. Tuning the turbulent viscosity pa- 
rameter will best be done by fitting to amplitudes of models with 
the best-known masses such as, e.g., those Idnematically associated 
with the thick disk. 
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